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Abstract 

Using T(2S') and T{3S) data collected with the CLEO III detector we have searched for decays of 
Xbj to final states with open charm. We fully reconstruct D^ mesons with p^jo > 2.5 GeV/c in three 
decay modes (K~Tr~^, K~Tr~^7r^, and K^Tr^ir^-K^) in coincidence with radiative transition photons 
that tag the production of one of the Xbj{nP) states. Significant signals are obtained for the two 
J = 1 states. Recent non-relativistic QCD (NRQCD) calculations of Xbj{nP) -^ ccX depend on 
one non-perturbative parameter per Xbj triplet. The extrapolation from the observed D^X rate 
over a limited momentum range to a full ccX rate also depends on these same parameters. Using 
our data to fit for these parameters, we extract results which agree well with NRQCD predictions, 
confirming the expectation that charm production is largest for the J = 1 states. In particular, for 
J = 1, our results are consistent with ccg accounting for about one-quarter of all hadronic decays. 



I. INTRODUCTION 

The six known Xbj{nP) P-wave bound states of a bottom quark (6) and its antiparticle b 
are labeled by their total angular momentum J = 0, 1, 2 and radial quantum number n = 1, 2. 
Their decays provide a place to test predictions based on Quantum Chromodynamics (QCD), 
which describes the strong interaction between quarks in the Standard Model of particle 
physics. While strong coupling prevents QCD at low energies from being treated with naive 
perturbation theory, specialized calculational techniques have been developed and applied 
with general success. In the bb system of states, one can study both transitions among the 
various quantum states, which also include the S'-wave T states, or else study decays which 
are initiated by annihilation of the quark- ant iquark pair. Although the Xbj states have been 
known for many years and there have been several studies of their transitions to other bound 
states in the bb system, there are no published annihilation decay branching fractions. This 
Article reports the first observation of some of the inclusive decays of the Xbj{^P, 2P) to D'^ 
mesons. 

In practice, one studies Xbj produced via the radiative transitions T{mS) -^ IXbjijiP) 
from T mesons produced directly at e^e~ colliders. The transition photons are typically 
used to tag XbJ events. Most of the XbJ radiative decays to the T states are well-measured 
[1]; the largest branching fraction is quite substantial, about 35%. Small, 0(1%), hadronic 
transitions to other bottomonium states, Xbi,2{'^P) — ^ 7r7rxfei,2(lP) and X6i(2P) — > OJ T(15'), 
have recently been observed [2] . The remainder of the decays are expected to be dominated 
by bb annihilation. Positive C-parity forbids decays via a single photon; the leading process 
is annihilation into two gluons. For the J = 1 state, decay into two on-shell gluons is 
forbidden [3]; instead, this state decays preferentially via qqg. While the J = 0,2 decay 
widths are dominated by this gg process, they also have a small admixture of qqg. 

We observe bb annihilation as a decay into lighter hadrons and are seeking to determine 
whether production of charm hadrons is suppressed or not. It is well-known that in contin- 
uum hadronization {e^e~ ^ 7 — > qq) that charm is not suppressed, while in ggg decays of 
the T(l^), an upper limit on £)*+ production of i3(T(15) -^ ggg -^ D*+X) < 1.9% (90% 
CL) indicates significant suppression [4]. 

The earliest calculations of inclusive charm (ccX) production from bottomonia focused 
on T — > ggg decays, giving estimates of a few percent [5]. It was soon pointed out that 
while production of ccX is predicted to be suppressed in gg hadronization, it is not expected 
to be suppressed in qqg hadronization [6]. Since the gg process is absent for the XbiinP) 
states, they should have higher branching fractions to ccX. These first calculations exhibited 
infrared divergences manifested as logarithms of the binding energy which were estimated 
in terms of a confinement radius. The predicted ratios of branching fractions are [6] Rj = 
^{Xbj -^ gg,qqg -^ ccX)/B{xbj -^ 99, qqg) = 6%, 25%, and 12% for the J = 0, 1, and 2 
states, respectively. The predictions were independent of the radial quantum number, n. 
The 25% branching fraction for J = 1 corresponds to equal rates for all accessible quark 
flavors q in qqg. 

With the development of non-relativistic QCD (NRQCD) techniques [7], a proper treat- 
ment of the infrared divergences was given and thus much improved calculations became 
possible. However, initial work [8] on bottomonium decays approximated final-state quarks 
as massless. Recently, this was remedied, and detailed NRQCD calculations of massive 
charm production in Xbj decay have been performed [9] . Decay rates are expressed in terms 
of one non-perturbative parameter per Xbj triplet: ps = ml{0^) j {Oi) where Oi (Cg) is a 



particular color-singlet (color-octet) four-quark operator [8, 9] and mf, is the one-loop pole 
mass, nib — 4.6 GeV/c^. All of the n-dependence in these calculations is contained in pg, 

(c) 

and Rj is found to increase monotonically with increasing pg. For illustrative purposes, 
we choose a common nominal value of pg = 0.10, which gives Rj = 5%, 23%, and 8% for 
the J = 0, l,and 2 states, respectively. These results are in general agreement with the 
older calculation cited above. In particular, charm production is expected to be largest for 
the J = 1 states. Not only the predicted Rj , but also the efficiency of our applied D^ 
momentum cut, depend on pg. We thus fit for pg in the context of the NRQCD results in 
order to interpret the consistency of our results with theory. 

To summarize, we observe charm production by observing D^ mesons in Xbj decays. We 
thereby hope to test predictions for the branching fractions, especially the expectation that 
the largest branching fractions will come from the J = 1 states due to the dominance of 
qqg decays when gg is absent. Sections II- VII present our experimental results for inclusive 
decays of XbJ to D^X, with a D^ momentum cut. Section VIII makes the connection 
between these measurements and the theoretically-predicted total rate of ccX production, 
Rj ■ Section IX summarizes our conclusions. 

II. THE CLEO III EXPERIMENT AND DATA SETS 

We use data collected with the CLEO III detector [10] at the Cornell Electron Storage 
Ring (CESR). Charged particle tracking is provided by a four-layer silicon tracker and a 
47-layer drift chamber [11] covering 93% of the solid angle. Particle identification (PID) is 
performed via specific ionization measurements (dE/dx) in the drift chamber supplemented 
by a Ring-Imaging Cherenkov detector (RICH) [12] which covers 80% of the solid angle. 
Photons are detected using an electromagnetic calorimeter consisting of 7784 CsI(Tl) crystals 
[13]. All of these detector elements are immersed in a 1.5 T solenoidal magnetic field. 

We use CLEO III data samples of 0.65, 1.27, and 1.40 fb^^ at the T{1S),T{2S), and 
T(3S') resonances, corresponding to 13.0, 9.4, and 6.1 million T mesons produced, respec- 
tively. In addition, data were also collected about 25 MeV below each resonance: we analyze 
0.14, 0.43, and 0.16 fb~^ from below the T(IS'), T(2S'), and T(3S') resonances, respectively. 
We do not use a direct off-resonance subtraction, but rather use these samples to constrain 
background shapes. 



III. EXPERIMENTAL TECHNIQUE 

This analysis includes all six known Xbj{nP) states: J = 0, 1, and 2 and n = 1 and 2. The 
XbJ states produced in radiative T decays are tagged by transition photons from T -^ JXbJ 
decays; the XbJ yields are obtained from fits to E^ spectra. We then fit E^ spectra from 
events with a D° candidate in the signal mass region, using D^ mass sidebands to remove 
combinatorial background under the D^ signal peak. After correcting for D° efficiencies 
and branching fractions, the ratio of these two inclusive yields determines the fraction of 
XbJ decays with a true D^ (above our D^ minimum momentum requirement). The photon 
efficiencies, numbers of initial T(nS'), and many associated systematic uncertainties largely 
cancel. 

We finally apply some small corrections to obtain the rate for direct production of D^ 
mesons in Xbj decays. Direct denotes the exclusion of charm production in decays of other 



bottomonium states produced by tranistions from our initial XbJ (for example, via 'y,7nr,u! 
transitions). Our focus is on direct D^ production via hadronization of xtj -^ 99, (1(19 decays 
only, and not on transitions to other bb states which subsequently decay to D^X. 

IV. EVENT SELECTION 

We first select events with transition photon candidates with energies between 3.50 < 
ln(£'^ [MeV]) < 5.70 (33 < Ey < 299 MeV). Only showers in the barrel calorimeter, 
I cos^l < 0.8, that are isolated from charged tracks are considered. Hadronic shower frag- 
ments are suppressed by vetoing any candidate photon shower that has a charged track 
pointing anywhere in the candidate's "connected region": this is a contiguous group of ad- 
jacent crystals with the energy deposition in each crystal, E^tai, satisfying E^tai > 10 MeV. 
An additional requirement on the fraction of energy deposited in the central 3x3 square of 
a 5 X 5 square, E9/E25, is applied. We use an energy-dependent E9/E25 criterion to select 
soft transition photon candidates, while photons later used in forming 7r° candidates, both 
as a veto and as D^ decay daughters, must satisfy the requirement of E9/E25 > 0.85. 

Photon background in the T — > 'JXbJ transitions is dominated by n^ decay products. To 
suppress this background, we reject photon candidates that, when combined with any other 
photon, form a vr*^ candidate that has an invariant mass within three standard deviations 
of the nominal vr" mass and a lab-frame opening angle between the two photons satisfying 
I cos^^^l > 0.7. 

For D^ reconstruction, we select well-measured tracks consistent with originating from 
the interaction point. These tracks must have an impact parameter of less than 5 cm with 
respect to the interaction point along the beam direction, and less than 5 mm with respect 
to it in the transverse plane. Charge-conjugate final states, D^X, are also included and 
are implied in the remainder of the paper. Candidate D^ mesons are reconstructed via 
three decay modes: K'n^, K~7i~^7i^, and i^^7r~7r'^7r+. For charged pion and kaon selection, 
particle identification combines RICH measurements with dE/dx in a momentum-dependent 
manner. The dE/dx information is expressed as cr^^, the number of standard deviations 
between measured and expected ionization for the vr, K hypothesis. The track-dependent 
dE/dx resolution used to normalize (Tt^^k includes dependencies on velocity, cos^, and the 
number of hits used for dE/dx. RICH information is characterized with a likelihood L; we 
use Ctt^k as shorthand for —2ln Lj^^k- When used, the RICH information is combined with 
dE/dx into one combined separation variable as: Ax'i k = ^■n,K — ^k,tt + {.(^t\Y ~ ((^IFttY- 
The first (second) subscript is chosen for tt {K) identification. We also impose requirements 
on the number of detected Cherenkov photons, n'Z'^ , for either the tt or K hypothesis in the 
RICH detector. 

Momentum dependence in the use of the RICH is motivated by the Cherenkov threshold 
for kaons and the need for tracks to have sufficient transverse momentum to reach the RICH 
detector given their curvature in the magnetic field. All pion candidates must satisfy |o"^^| < 
3. Pion candidates with p < 0.50 GeV/c are accepted with that criteria alone, but additional 
requirements are added for some higher-momentum candidates. If 0.50 < p < 0.65 GeV/c 
and n^ > 2, we also require Ax'i < 0. Candidates with p > 0.65 GeV/c must satisfy both 
n^ > 2 and Axl < 0. 

Kaons are identified in an analogous manner to pions, with three additional criteria. 
First, kaon candidates must satisfy p > 0.18 GeV/c. Kaons lose more energy in the inner 
detector than pions, and tightly curling tracks are poorly reconstructed. Second, if the 
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FIG. 1: Sum of A'^vr^, A^vr^vr*', and K~'k~'k^'k^ invariant mass distributions obtained for T(2S') 
(a) and T(3S') (b) data. The shaded areas correspond to the signal region and the two background 
side-band regions defined in the text. 



K 7r^7r° are reconstructed from pairs of photons 



track momentum is greater than 0.60 GeV/c, then the track must also be within the RICH 
fiducial region, | cos^| < 0.80; this ensures good rejection of the more numerous pions as the 
dE/dx separation degrades. Finally, when RICH information is available, a tighter criterion, 
Axl^ < —10, is used compared to that employed for pions due to the relative abundance of 
pions over kaons. 

The 7r° meson candidates from D^ 
with an invariant mass within 2.5 standard deviations of the nominal vr" mass. These 
candidates are then kinematically constrained to the 7r° mass. For the K^tx^-k^ mode, the 
precision is improved with an additional requirement on the candidate's location in the 
Dalitz plot. Our criteria retains the 70% of decays from the most densely-populated regions 
of phase space (based on previous measurements [14]). 

In order to avoid the large combinatorial backgrounds under the D^ signal at lower 
momenta, only candidate D^ momenta p^o > 2.5 GeV/c are accepted. Figure 1 shows the 
sum of the K~n~^, K^ti^tt^, and K'Ti^-n-^Tr^ invariant mass distributions, Knii {n = 1, 2, 3), 
obtained from T{2S) and T(35') data for events also containing transition photon candidates. 
The D^ signal region is defined as the K~n~^, K~7i~^7[^, and K^Tf^Tr^ir^ invariant mass 
interval ±2.5(Jm (using a mode-averaged am — 0.0075 GeV/c^) from the nominal D^ mass, 
m£,o [1]. The D° "sideband" regions, each with a width of 2.5am, are located symmetrically, 
between 7.5 am and 10.0 am on either side of the nominal D^ mass. 



V. FITS TO THE PHOTON ENERGY SPECTRA 

We first measure the total number of Xbj tagged with an observed transition photon by 
fitting the inclusive E^ spectrum. Photon peaks from inclusive T(2S') -^ 7Xbj(l-P) and 
T{3S) -^ ■yxbji'^P) transitions are evident in Fig. 2. 
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FIG. 2: Fits to the T{2S) (left) and T(35) (right) inclusive photon energy spectra. The data are 
shown as dots; the fits are shown as the histograms; the dashed lines represents the total fitted 
background. Note the suppressed zero on the vertical axis. Nominal photon peak locations for 
transitions to the Xbj{^P) (on the left) are 111, 130, 164 MeV/c^ (for J = 2, 1,0, respectively) and 
for transitions to the Xbji'^P) (on the right) are 87, 100, 123 MeV/c^ (for J = 2, 1, 0, respectively). 



We use T(15') resonance and T{nS) off-resonance data to model the photon background 
in the E^ spectra [15]. The off-resonance data are observed to have indistinguishable spectra 
in our energy region and thus the three samples are combined to increase statistics. The 
T(IS') on-resonance and T{nS) off-resonance shapes are also quite similar, and we initially 
fit with two independent normalizations to peak-free regions of the photon energy spectrum. 
The regions are defined by 3.50 < \n{E^ [MeV]) < 3.70 (33 MeV < E^ < 40 MeV) and 5.55 < 
ln{E^ [MeV]) < 5.70 (257 MeV < E^ < 299 MeV) and the fit results are used to then fix 
the relative normalization of these on- and off-resonance samples for subsequent signal fits. 

When fitting the full photon energy spectra to extract signal yields, only one overall 
normalization parameter for the background is varied. We find, however, that the fit quality 
is acceptable only after the inclusion of first- (1-P) or second-order (2P) polynomials to allow 
small smooth adjustments of the background shape. The fit also includes signal contributions 
from the three dominant El transitions, T(2S') -^ ■jXbji^P) or T{3S) -^ 7Xfej(2P), as 
appropriate. The X6j(lP) and xtji'^P) signal peaks are described by a so-called Crystal Ball 
line shape [16] with fixed asymmetry parameters, a and n. This line shape is a Gaussian, 
described by a peak energy Ep and resolution aE, matched with the constant c onto an 
asymmetric low energy tail, l/{Ep — E + c)", at an energy Ep — aae- We obtain Ep from 
published results [1] and use the values a = 0.84 and n = 25.8. The values of sigmaE/E 
depend on E, varying from 5.4% to 3.9% as the energy of the six transition lines increases. 
This E dependence is determined from Monte-Carlo studies, but the overall scale of the 
resolution is adjusted based on fits to data. In addition to the dominant T(35') -^ 7X6j(2P) 
transitions, the fit to the T{3S) spectrum includes the lines due to Xbji.'^P) —>■ 7T(2S') 
cascades. The fit results are displayed with the data in Fig. 2 and tabulated in Tables I and 
II. 

Photon energy spectra for events with D° mesons are obtained by subtracting the ln(£'^ 
[MeV]) spectra associated with the Knn {n = 1,2,3) D° sidebands from the D^ signal 
region. The ln{E.y [MeV]) distributions and the fits for the T(2S') and T{3S) data are 
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FIG. 3: Fits to the T(25) (left) and T(35) (right) photon energy spectrum obtained for events 
with D^ mesons. The data are shown as dots; the fits are shown as histograms; the dashed hnes 
represents the total fitted background. 
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FIG. 4: Energy spectrum for background-subtracted T(25) -^ ^Xbji^P) (left) and T{3S) -^ 
lXbj{'^P) (left) photon lines obtained for events with D^ mesons. The data are shown as dots; the 
fit is shown as the solid line. Individual contributions from the signal T{mS) -^ "fXbj{nP) lines 
are shown as dashed-line peaks. 

presented in Fig. 3. The J = 1 lines are the most pronounced. Photon background shapes 
for these spectra are the same as for the T{2S) and T(3S') inclusive photon analysis, except 
that an acceptable fit quality is obtained without the addition of low-order polynomials, and 
they are omitted. The background-subtracted photon spectra are presented in Fig. 4 and 
fit results are tabulated in Tables I and II. 



VI. MEASUREMENT OF XbJ -^ D°X (po" > 2.5 GeV/c) RATES 



K TT^TT^^K 71 11 



The yields of events with xtj and Z)° mesons (-D° ^ K tt^ , A vr ' tt", A vr tt ' vr ' j m- 
clude non-direct Xbj decays which must be subtracted. Non-direct Xbji^P) decays to D^X 



include T(25') -^ 7Xw(l-P); Xbj{^P) -^ I'^i'^S) decays where D'^ mesons are then produced 
in T(IS') annihilation into ggg, gg'~f, and 7. 

Non-direct Xbji'^P) decays to D^X similarly include production of bottomonium states 
which in turn may decay to D'^X. Known processes include T(IS') produced via T(3S') — > 
'jXbji'^P) followed by 

. xw(2P) ^7T(25);T(25) ^ (nn , n^ , r])T (IS) 
. xw(2P) -^ 7T(25); T(25) -. iXbAlP); XbAlP) ^ lT{lS) 
• Xbj{2P) ^ nnxui^Py, Xbj{lP) ^ iTilS) 
and Xbj{^P) produced via T(3S') -^ 7X6j(2P) followed by 
. Xbj{2P) ^ TrnxuilP) 
. XbA2P) -. 7T(25); T(25) ^ 7Xw(l^) 

and T(2S) from T(3S) -. 7Xw(2P); Xw(2P) ^ 7T(25). 

Yields for events with D^ mesons from direct Xbj{^P) decays are calculated by correcting 
raw yields from the T{2S) data with a non-direct rate determined using known branching 
fractions [1] and an T(IS') -^ {999,991,1) -^ D^X rate for p^jo > 2.5 GeV/c of 2.60 ± 
0.50% [17]. We estimate the numbers of these non-direct events as 16±9, 191±58, and 125± 
34 for J = 0, 1, and 2, respectively. Corresponding estimates of the non-direct backgrounds 
for Xbj{2P) -^ D^X in the T(3^) data are 53 ± 24, 392 ± 70, and 311 ± 50 for J = 0, 1, and 
2, respectively. We account for the fact that prompt production of D^X from T(25') differs 
from that from T(IS') due to the different mixture of decays mediated by ggg,ggi, and 7. 

Yields for inclusive Xbj production, total Xbj with D^ mesons, and Xbj with directly- 
produced D^ mesons, are summarized in Tables I and II. In addition, we list a correction 
due to a small observed curvature in the Kmr mass spectra leading to a small residual 
background of true photons and fake D^ mesons, since our sideband subtraction assumes a 
fiat background. 

The direct XbJ yields, A'^^^ '*'', from A^t initial T produced are: 



Kf^ = ^r e, B{T ^ 7Xw) Bixbj - gg, qqg ^D'X)Y: e.i3.(D°), 

where e^ is the 7 detection efficiency and the last factor I]eji3j(D°) is a sum over the three 
Knii decay modes of the D^. The observed number of inclusive Xbj decays is given by 

N'^) = Nre,B{T^^Xbj). 

Our main results, the branching fractions B{xbj —^ gg, QQg -^ D^^), are obtained from the 
two previous equations as: 



^iXbJ -> 99, 119 -^ D^X) 



XbJ 



NllfEe^B.iD'^y 



where the photon efficiency, e^, and sample size, iVx, both cancel. For determination of the 
D^ detection efficiencies, Monte-Carlo simulation of continuum cc events (based on Jetset 

9 



7 [18]) were used, since this sample is expected to approximate the jet-hke events from 
the XbJ -^ ccg decays. We find that the efficiency is consistent with being independent of 
momentum in the p£,o > 2.5 GeV/c range. 

Based on detailed comparisons of particle identification in our data and Monte-Carlo 
simulations, we conclude that small efficiency corrections are needed. The Knir modes 
receive adjustments of fxf^, where /k = 0.95(0.99) and U = 0.99(1.01) for T{2S) (T(35)) 
data. The Xbj -^ D^X decay rates for p^o > 2.5 GeV/c are presented in Tables I and II. 

TABLE I: T(2S) -^ 7Xbj(l-f) i-J = 0;1)2) transition yields and Xb ~^ 99iQQ9 ~^ D^X rates, for 
P£)() > 2.5 GeV/c. Errors shown are statistical only. 

Final state Xfeo(l^) Xbi(lP) Xb2{lP) 

Mfi 166860 ±5988 363825 ± 6793 379457 ± 7243 

XbJ 

N^l] (raw) 501 ± 303 2561 ± 346 1207 ± 360 

D^ sideband correction llib5 60 ±6 57 ±7 

non-direct Z)° 16 ± 9 191 ± 58 125 ± 34 

^Xb/'"" (direct) 474 ± 303 2310 ± 351 1025 ± 362 

B{xbjO-P) -^ 99,qqg -^ D^X) 5.63 ± 3.61% 12.59 ± 1.94% 5.36 ± 1.90% 



TABLE II: T(35') -^ jXbji'^P) {J = 0, 1,2) transition yields and Xb ~^ 99^119 ~^ D^X rates, for 
Pdo > 2.5 GeV/c. Errors shown are statistical only. 

Final state Xbo{2P) Xbi{2P) Xb2{2P) 

Mfi 219773 ±5201 491818 ± 5197 524549 ± 5628 

XbJ 

N^°j (raw) 565 ± 341 2757 ± 366 477 ± 370 

L»° sideband correction 39 ± 7 122 ± 7 122 ± 7 

non-direct Z)0 53 ± 24 392 ± 70 311 ±50 

N^Uir (direct) 473 ± 342 2243 ± 373 44 ± 373 

B{xbj{2P) -^ gg,qqg -^ D^X) 4.13 ± 3.00% 8.75 ± 1.47% 0.16 ± 1.37% 



VII. SYSTEMATIC UNCERTAINTIES ON THE BRANCHING FRACTIONS 

Systematic uncertainties on the six measured branching fractions are primarily of two 
types. The first are uncertainties in D° reconstruction; these affect each of the six Xbj states 
equally and are summarized in Table III. The next are uncertainties related to our photon 
yields, both in terms of efficiencies and yield extractions. These often differ for the six Xbj 
states, and are summarized in Table IV. In the remainder of this section we detail the 
sources of the uncertainty estimates presented in the aforementioned Tables. 
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The first three entries of Table III involve efficiencies for track-finding, 7r° reconstruction, 
and particle identification algorithms. Since the composition of the three D^ final states 
differ, we take a linear weighting of the uncertainties across D^ modes. The weights used 
are Wi = eiBi/ Ej^jEj, yielding 0.25, 0.34, and 0.41 for L)° -^ K-n+, D° -^ R-n+n^, and 
D^ -^ A'^TT^TT+TT"'", respectively. 

The systematic uncertainty in track-finding is obtained by studies of the difference be- 
tween data and Monte-Carlo simulation. We assign a 1.5% uncertainty per track, which 
gives a net uncertainty of 4.2% after weighting across D^ decay modes. 

We assess the uncertainty in vr^-finding at 5% per tt". Taking into account the weight of 
the D^ -^ K'71^71^ mode, the net Tr^-finding systematic uncertainty is 1.7%. 

Systematic uncertainties in kaon and pion identification are obtained by comparing data 
and Monte-Carlo efficiencies. We obtain 2% (1%) uncertainties per K (vr) which yield a net 
4.0% systematic uncertainty, averaged over D^ modes. 

The systematic uncertainty on the D^ -^ i^~7r+7r° efficiency due to selection on the Dalitz 
region is obtained by comparing the inclusive yield changes in data compared to Monte-Carlo 
simulations as the selection efficiency is varied. As a result of this study, and accounting for 
the fraction of D^ candidates found via this decay mode, we assign 1.0% as our total Dalitz 
region selection uncertainty. 

For evaluation of systematic uncertainties related to the D^ momentum requirement, the 
Poo requirement was varied. Events were selected for three values of the D^ momentum re- 
quirement (> 2.2, > 2.5, and > 2.8 GeV/c). We assign a 1.7% branching fraction uncertainty 
due to this source. 

To study possible effects of the event shape and environment on the D° detection effi- 
ciency, different models of signal Monte-Carlo and continuum Monte-Carlo events are an- 
alyzed. Results indicate a 3.0% uncertainty of the efficiency for the event-shape changes 
explored. 

Systematic uncertainties related to the definition of the D^ signal and sideband regions 
are obtained by varying the corresponding mass windows. This also includes uncertainty 
due to a nonlinear background shape under the D° signal. The total systematic uncertainty 
is determined to be 2.5%. 

The total uncertainty in the D° efficiency is 7.5% for each Xbj state, as noted in Table 
111. We now turn to the photon-related systematic uncertainties presented in Table IV. 

TABLE III: Relative systematic uncertainties on measured branching fractions from sources affect- 
ing the D^ efficiency. 



Source 




Uncertainty (%) 


Tracking: 1.5%/track 




4.2 


7r° efficiency: 5%/7r° 




1.7 


PID: 2%/K^, l%/7r± 




4.0 


Kmr^ Dalitz requirement 




1.0 


Momentum dependence 




1.7 


Decay model effects on De- 


' efficiency 


3.0 


selection of events with a 


D^ 


2.5 


Total D^-related systematic uncertainty 


7.5 
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To verify that the photon efficiency largely cancels in our analysis, the difference of photon 
efficiencies between inclusive events and those with a -D° candidate is studied using Monte- 
Carlo samples. We find that the relative photon efficiency difference between spherical ggg 
events and jet-like qq events is about 6%. In our case, we are concerned about the difference 
between generic Xbj events and those having a reconstructed D^ . Presumably the effect of 
this bias is smaller than that of the rather large overall event shape change between these 
two Monte-Carlo samples. We thus take 1/3 of the variation and assign a 2% uncertainty 
for all six Xbj states. 

For estimation of line-shape ffiting uncertainties we change the Crystal Ball line-shape 
parameters a and n by ±10% from their nominal values. This range is chosen as appropriate 
based on changes in ffi quality. We take the resulting branching fraction variations as 
systematic uncertainties, ranging from 0.1% to 0.6%. 

The nominal ffiting ranges for photon energy distributions are 3.8 < hi{E^ [MeV]) < 5.5 
for T(2S') and 3.8 < hi{E^ [MeV]) < 5.7 for T(3S'). We vary the lower and upper limits of 
the ffiting regions from 3.50 to 3.70 and from 5.50 to 5.70 . Variations in our results suggest 
uncertainties from 0.3% to 0.6%. 

As mentioned above, the photon background shape consists of two components: the reso- 
nant and off-resonance photon spectra used to estimate the background shapes in the T(2S') 
and T(3S') photon energy distributions. We varied scaling factors for the photon back- 
ground components and changed the T(IS') resonance and the T(2S'), T(35') off-resonance 
contributions in the photon background shape. Also, in the ffi of the T(2S') and T(3S') 
inclusive photon energy distributions, we used additional background components to obtain 
a better ffi quality. First, second, and third order polynomials are tried as extra components 
in addition to the T(15') on-resonance and the T(2S') and T(3S') off-resonance background 
shapes. We estimate systematic uncertainties due to such choices at levels ranging from 
0.5% to 1.6%. 

Our nominal ffi uses logarithmic binning of energy ln(i5^ [MeV]). We changed the log- 
arithmic energy scale to linear binning, with 1 MeV energy bins. The photon background 
shape was left unchanged. We assign from 0.2% to 1.7% uncertainties on our branching 
fractions based on the stability of our results. 

The T (35") photon energy spectrum includes T(2S') -^ jXbj{^P) transition lines at similar 
energies. To estimate systematic uncertainties on the B{xbj{'^P) -^ -D^AT), we include these 
lines in the ffi to the T(3S') inclusive photon spectrum and the photon spectrum for events 
with D^ mesons. Estimated systematic uncertainties varied from 0.2% to 1.5%. 

In Table IV, we summarize the systematic uncertainties associated with 7 detection and 
ffiting for each of the six Xbj lines. Note that these uncertainties apply to the raw yields, 
before any subtractions are made. 

We also performed several simple cross-checks to investigate the stability and consistency 
of our results. These included splitting the datasets into two subsets, varying selection 
criteria, and comparing yields in individual D^ decay modes. All of these tests produced 
consistent results. 

Our final results for p^o > 2.5 GeV/c are given in Table V. Upper limits are given for 
modes without significant signals, but central values for those modes will be needed for ffis 
later. 



12 



TABLE IV: Relative systematic uncertainties on measured branching fractions due to sources 



related to the E^ distributions. 









Uncertainty (%) 






Source 


XboilP) 


XbiilP) 


Xb2{lP) 


Xbo{2P) 


Xbi{2P) 


Xb2{2P) 


7 efficiency cancellation 


2.0 


2.0 


2.0 


2.0 


2.0 


2.0 


Line-shape fitting 


0.6 


0.1 


0.4 


0.5 


0.1 


0.5 


Fitting range 


0.5 


0.3 


0.4 


0.6 


0.3 


0.5 


Background shape 


1.4 


0.6 


0.9 


1.6 


0.5 


0.9 


7 energy binning 


1.4 


0.2 


0.5 


1.7 


0.3 


0.6 


T{2S) ^ -fXbAlP) lines 


- 


- 


- 


1.5 


0.3 


0.2 


Total 7 systematic uncertainty 


2.9 


2.1 


2.3 


3.5 


2.1 


2.4 



TABLE V: Summary of measured branching fractions (or upper limits) for B{xbj{nP) -^ gg, qqg -^ 
D^X) with the requirement that p^o > 2.5 GeV/c. The uncertainties are statistical and systematic, 
respectively. 



State 


B{xb. 


i{nP) ^ gg,qqg^ D' 


^X) (%) 90% CL UL (%) 


XboilP) 




5.6 ±3.6 ±0.5 


< 10.4 


XbiilP) 




12.6 ±1.9 ±1.1 




Xb2{lP) 




5.4 ±1.9 ±0.5 


<7.9 


Xbo{2P) 




4.1 ±3.0 ±0.4 


<8.2 


Xbi{2P) 




8.8 ±1.5 ±0.8 




Xb2{2P) 




0.2 ±1.4 ±0.1 


<2.4 



VIII. INTERPRETATION 

We observe significant production of D^ mesons from both the Xbii^P) and Xfei(2-P) 
states. There is evidence of a signal for Xb2{^P), while data for the other three states are 
inconclusive. For each triplet, we observe the largest branching fraction for the J = 1 states, 
as expected. 

The NRQCD calculation mentioned earlier [9] makes predictions for the total ccX pro- 
duction rate, Rj, as a function of one non-perturbative parameter, pg, per Xbj triplet. We 
would like to convert our measurement of the inclusive D^X rate, with a minimum momen- 

(c) 

tum requirement, into an experimental value for Rj . However, this conversion also depends 
on ps, since this parameter affects the momentum spectrum of the D^ mesons and hence the 
efficiency of our minimum momentum requirement. We use six branching fraction results to 
determine two best-fit values of pg (one per triplet). Our experimental results for Rj are 
based on these best-fit values and clearly depend on our use of the NRQCD calculation. 
We first discuss the details of how to relate our measurements to the inclusive ccX rate 
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and then present our extraction of the pg parameter and experimental values of Rj ■ Three 
factors will combine to cause our extracted Rj to be larger than the directly-measured 
branching fractions in Table V. We only see some of the D^ spectrum, not all charm 
appears as D'^, and Rj is normalized to the number of Xbj that decay via annihilation, not 
the total number produced. Only one factor works in the other direction: Rj measures 
ccX production, and either charm quark may form a D^ . 

Suppressing the Xtji^P) radial quantum numbers for simplicity, we have: 



R 



J 



^{XbJ -^ 99, qq9 -^ ccX) B{xbj -^ 99, 1(19 ^ D^^, Pd^ > 2.5 GeV/c) 



^iXbJ -^ 99, m9) /2.5 /do B{xbj -^ 99, QQ9) 



where the right-hand side contains our directly-measured branching fraction with three ad- 
ditional factors which we now explain. 

First, we must divide by B{x ~^ 99, QQ9) such that the final branching fraction is normal- 
ized to only gg, qqg decays of the Xbj since this is the normalization used for the theoretical 
prediction. These branching fractions are calculated as 1 — I]j iS^, where the sum extends 
over all known transitions of a given XbJ to other bottomonium states [1] . 

Next, we divide by /2.5, the fraction of the D^ spectrum expected to be above our 
2.5 GeV/c D^ momentum requirement. This is obtained from the results of Ref. [9], and it 
depends on the value of pg and knowledge of the charm fragmentation function [19]. 

Finally, we must divide by the number of D^ mesons expected per ccX event: /^o = 
1.11 ± 0.08. This number is itself the product of four factors. The first is a factor of two 
to account for the two quarks, each of which may form a D^ . The next two factors account 
for all seven weakly-decaying C = 1 states D°, D^, Ds, Ac, S^, S°, and Q^, relative to 
the measured D^ yields. The fraction of D^ compared to the total of D° -|- D~^ + Dg + 
Ac, N{D'^)/[N{D^) + N{D+) + N{Ds) + iV(Ac)] = 0.574 ± 0.041, is obtained from e+e" 
fragmentation data [19]. An additional factor 0.98 ± 0.01 then accounts for the omitted 
H+, S°, and fi° states. This is estimated from the A^ fraction of N{Ac)/[N{D°) + N{D+) + 
N{Ds)] = (8.1 ±2.1)% in [19] (with an added uncertainty from knowledge of B{Ac — ^ pKtt)), 
combined with a theoretical suppression of order 10% due to the additional strange quark 
popping needed to form the omitted states. The fourth factor of 0.99 ± 0.01 accounts for 
charmonium states, which here include those states below open-flavor threshold at ^/s = 
2Mbo: J/ij,ij{2S),r]c,r]c{2P),Xcj,hc. We estimate iV(openc)/[A^(openc) + 2A^(cc)] ~ 1 - 
2A^(cc)/A^(openc) ~ 1 — 2B{ccX — * charmonia) based on the production rate of J/ip in 
e~^e~ fragmentation [20] and the branching fractions to charmonium in T(IS') decays [1]; 
these processes show that charmonium is rare in both 7 and ggg hadronization. We are not 
sensitive to errors at the 1% level and choose a conservative uncertainty to accommodate 
unmeasured charmonium states. The various factors required for the six XbJ states are 
summarized in Table VI. 

With these factors in hand, we fit our data for the D^X branching fractions with p£,o > 
2.5 GeV/c to the NRQCD predictions [9] and extract pg, the ratio of color-octet to color- 

(c) 

singlet matrix elements, in Xbj decays. Recall that both /2.5 and Rj depend on pg and 
that /2.5 depends on fragmentation functions. For each value of pg, we may convert our 
directly measured branching fractions into extracted values for Rj in the context of this 
NRQCD calculation (which includes the assumption that e~^e~ charm fragmentation data 
is representative of our charm fragmentation). The best value of pg is obtained from a fit 

(c) 

which finds the best agreement between the predicted and extracted Rj ■ 
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TABLE VI: Summary of factors used to relate our measured D^X branching fractions to Rj, 
which measures the total ccX rate. The values of /2.5 are evaluated at the independently fitted 
best values of pg for each triplet. 

Factor XbojlP) XbijlP) Xb2{lP) Xbo{2P) Xbi{2P) Xb2{2P) 

Six -^ 99, qq9) 0.97 ± 0.03 0.65 ± 0.08 0.78 ± 0.04 0.93 ± 0.07 0.68 ± 0.04 0.75 ± 0.03 
/2.5 0.54 0.70 0.63 0.45 0.46 0.47 

fjjo 1.11 ± 0.08 1.11 ± 0.08 1.11 ± 0.08 1.11 ± 0.08 1.11 ± 0.08 1.11 ± 0.08 

1/(/d"/2.56) 1.70 ± 0.13 1.97 ± 0.28 1.83 ± 0.16 2.15 ± 0.23 2.89 ± 0.28 2.56 ± 0.21 



We fit separate ps values for each triplet by minimizing a x^ which has one term for 
each of the three states. Each term in the x^ is formed from the square of the deviation 

(c) 

of the predicted and extracted Rj values, normalized by the errors on the extracted value. 
Note that both the predicted and extracted Rj values depend on pg. Correlated systematic 
uncertainties on the branching fractions are incorporated into the covariance matrix used 
to evaluate the x^ ^^ o^^ fits. We find, however, that results are insensitive to correlations 
due to the dominance of statistical errors. The best-fit values are ps(lP) = O.lGOto^oIy and 
P8(2P) = 0.074:t'^om8 with x^(lP) = 0.40 and x^(2P) = 4.71, respectively, for 3 - 1 degrees 
of freedom each. The errors are larger for the IP states primarily due to the non-linear 
dependence of the branching fractions on pg: for larger pg, the branching fractions are less 
sensitive to changes in its value. 

It has been argued [21] that pg should be largely independent of radial quantum number. 
While we prefer not to assume such an equality, a joint fit to our branching fractions for 
both triplets obtains a best-fit common value of pg = O.OSGt^ois, with x"^ = 10.1 for 6 — 1 
degrees of freedom. 

Table VII lists the best-fit branching fractions, Rj , extracted from our data along with 
the best-fit NRQCD values, based on fits with separate pg parameters for each Xbj triplet. 
We also show the original 1979 calculations [6] for comparison. The third uncertainty is 
due to uncertainties in the branching fractions used to obtain B{x -^ gg, qqg) and the 
fragmentation data used to obtain f^o and /2.5. No systematic uncertainty is included for 
the accuracy of the theoretical calculations or the assumption that the e~^e~ fragmentation 
data is a valid model for our charm fragmentation since we do not know how to quantify 
such effects. Thus, while our primary results for the inclusive Xbj branching fractions into 
D^X with pf)0 > 2.5 GeV/c are model- independent, our results for Rj are clearly model- 
dependent. 



IX. CONCLUSION 

We report first measurements of the branching fractions for Xbj{^P,'2P) -^ D^X with 
Pd^ > 2.5 GeV/c. Our results are used to infer the total production of charm in Xbj decays, 
Rj in the context of a recent NRQCD calculation [9]. The results are in agreement with 
this calculation, as well as the older calculations [6]. Notably, our Rj values confirm that 
the largest branching fractions to charm correspond to the J = 1 Xbj states. 
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TABLE VII: Summary of extracted branching fractions (or upper limits) for Rj ■ NRQCD best-fit 
values use distinct p% values for each Xbj triplet. The original 1979 calculations [6] are also shown. 
The uncertainties are statistical, our systematic, and external systematic, respectively. 

State Rf (%) 90% CL UL (%) NRQCD Best Fit (%) Pred. from [6] (%) 

Xfeo(l^) 9.6 ±6.2 ±0.8 ±0.8 < 17.9 O 6 

Xfei(lP) 24.8 ±3.8 ±2.2 ±3.6 23.7 25 

Xb2(l^) 9.8 ±3.5 ±0.9 ±0.9 < 14.6 10.8 12 

X6o(2P) 8.7 ±6.4 ±0.9 ±0.7 < 17.7 4.9 6 

X6i(2P) 25.3 ±4.3 ±2.5 ±2.4 22.1 25 

X62(2P) 0.4±3.5±0.4±0.1 < 6.1 7.4 12 
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